Microscopic theory of thermal phase slips in clean narrow superconducting wires 
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We consider structure of a thermal phase-slip center for a simple microscopic model of a clean one- 
dimensional superconductors in which superconductivity occurs only within one conducting channel 
or several identical channels. Surprisingly, the Eilenberger equations describing the saddle-point 
configuration allow for exact analytical solution in the whole temperature and current range. This 
solution allows us to derive a closed expression for the free-energy barrier, which we use to compute 
its temperature and current dependences. 
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Recent advance in fabrication techniques brought ul- 
trathin superconductor wires into active experimental in- 
vestigation 

[3, 13, SI as they may play an important role 
in future electronics. These studies brought about the 
fundamental question of the limits till which thin wires 
retain a key property of superconductors, the ability to 
carry currents without dissipation. It was well estab- 
lished that in extremely narrow superconductors thermal 
fluctuations can occasionally drive a small part of the 
wire into a normal state. When there is only one chan- 
nel for a supercurrent to pass, these fluctuations inter- 
rupt its flow and cause a finite effective resistance. This 
phenomenon has been understood by Langer and Ambc- 
gaokar (LA) as thermally activated phase slips lo- 
cal abrupt switches of the system between its metastable 
states in which events phase jumps discontinuously by 
±27r, see also Ref. @ for a review. The corresponding re- 
sistance exhibits Arrhenius law-like temperature behav- 
ior, p tx exp(— SF/T), where SF is the energy barrier 
separating these states. In a phase-slip event the order 
parameter goes through the saddle-point configuration 
in which it changes sign as function of coordinate and 
vanishes at one point. The barrier for the phase slip, 
SF, is given by the energy of this saddle-point configura- 
tion. Later preexponential factor in resistivity has been 
calculated by McCumber and Halperin [j| and Langer- 
Ambegoakar-McCumber-Halperin (LAMH) theory has 
been confirmed experimentally on tin whiskers [9]. The 
dynamical theory of phase slips at large currents has been 
further extended in Refs. [ToL fill] . 

The mechanism of appearance of a finite resistivity in 
a thin wire can be described as follows [f\ . A finite volt- 
age drop across the wire means that the order param- 
eter phase difference </>i2 between the ends of the wire 
has to increase linearly with time. At first glance, this 
should inevitably lead to continuous increase of the su- 
percurrent flowing through the wire, which is obviously 
is inconsistent with steady state. The resolution of this 
apparent inconsistency is that the current relaxes via the 
phase-slip events, in which the order parameter vanishes 
at some point in the superconductor, allowing a phase 



difference between opposite sides of this point to change 
by 2tt. As a result of this event, a phase jump equal to 
2ir appears at one point of superconductor reducing the 
phase gradient in the rest part of the wire. To maintain 
a steady state, such events must occur with the average 
frequency 2e(V)/h. 

Phase slip events can be caused also by quantum tun- 
neling, which at sufficiently low temperatures can prevail 
over the thermal activation, and this possibility is being 
widely searched and debated [H, 0, S 0, 0] • Yet, the role 
and the very fact of observation of quantum phase slips 
remain a controversy. It seems now plausible that ther- 
mally activated phase slips may dominate down to very 
low temperatures and that quantum fluctuations do not 
add any significant corrections to the resistance in almost 
whole experimentally accessible temperature range. As 
original LAMH theory [g, |8| was developed within the 
framework of the Ginzburg-Landau (GL) approach, it 
holds, strictly speaking, only in the temperature interval 
near the critical temperature, T c , of the superconduc- 
tor transition. On the other hand, in sufficiently nar- 
row wires a finite thermally-activated resistivity can be 
observed down to temperatures which are significantly 
lower than T c , There is thus an experimental quest for ex- 
tension of the LAMH theory to low temperatures, where 
the phase slips have to be analyzed within the micro- 
scopic theory. 

In this Letter we investigate a simple model of a clean 
one-dimensional superconductor containing either only 
one conducting channel or a few degenerate channels. It 
turns out that this model allows for full analytical treat- 
ment of a free-energy saddle point configuration. An ex- 
act analytical solution of Eilenberger equation is found, 
both at zero and finite supercurrent. This solution allows 
us to derive a closed expression for the temperature and 
current dependent free-energy barrier. 

We consider the simplest case of a narrow clean super- 
conducting wire containing cither only one superconduct- 
ing channel or several degenerate channels, i.e. with the 
same Fermi velocity, vp ■ The wire diameter is assumed to 
be smaller than the London penetration depth; one can 
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thus neglect the magnetic field generated by the currents. 
We will find an order parameter configuration which cor- 
responds to free-energy saddle point in such a wire. As 
the saddle-point configuration is a stationary point in the 
parameter space, it is described by the same equations 
as the equilibrium state. In the quasiclassical approx- 
imation such superconducting wire is described by the 
gap parameter A(x) and the Green's functions, g, /, p 
depending on the coordinate along the wire x, the Mat- 
subara frequencies u n = 27rT(n+ 1/2) and the location 
on the Fermi surface, which is reduced to two points cor- 
responding to forward and back motion. The Green's 
functions are connected by the normalization condition 
g 2 + fp = 1 and obey the one-dimensional Eilenberger 
equations flij ]: 



dx 
df 

-v F -L -2u n f + 2A.g = 0, 
ox 

v F ^--2u n p + 2A*.g = 0, 
Ox 



(la) 
(lb) 
(1c) 



with v F = ±|i\f|. It is convenient to introduce new di- 
mensionless variables 

x -» x/i, A -» A/T c , T -> T/T c , 

where £ = v F /T c is the coherence length and T c is the 
critical temperature. In new variables the Eilenberger 
equations for forward direction, v F > 0, become 



dg_ 

dx 

dx 



dp_ 

dx 



hA/t-A*/ = 0, 
2u n f + 2Ag = 0, 
2uj n p + 2A*g = 0. 



(2a) 
(2b) 
(2c) 



These equations must be completed with the self- 
consistency equation, 



= KT^{f{u) n ,X,k F )) F , 



(3) 



and the formula for the supercurrent, 

3 = -2meN{0)T^2{v F g{u n ,x,k F )) F , (4) 

where N(0) = l/(irv F ) is the density of states. The 
summation is going over all Matsubara frequencies n = 
— 00...00. For the one-dimensional problem averaging 
over the Fermi surface, (. . .} F , reduces to summation over 
two points 

(/K, x, k F )) F = i (/(u> n , x, k F )+f(u n , x, ~k F )) , 
v F 

(v F g(uj n ,x,k F )) F = — (g(uj n ,x,k F )-g(uj n ,x,-k F )) . 



The Green's functions in the opposite directions are con- 
nected by relations f{—k F ) = [/^(Aif)] , JH^^f) = 
[f(k F )}* and g(-k F ) = [g(k F ))\ 

Let a finite current flow through the wire. Then the 
order parameter, A, and the anomalous Green func- 
tions / and p can be represented in the form of plane 
waves with slowly varying in space amplitudes, A, / oc 
gifcx. A* , p ex e~ lkx , where the wave vector k is de- 
termined by the supercurrent. Defining new variables 
f± = (/ ± p) /2, one rewrites Eqs. (!2aH2cp as 



dg_ 

dx 
df- 
dx 
df+ 
dx 



2iAjf + - 2A R f_ = 0, 

- 2A R g + 2JJ+ = 0, 

- 2iA ig + 2u;' n f- = 0, 



(5a) 
(5b) 
(5c) 



with An = ReA, Aj = ImA, and u' n = ui n + ik/2. Using 
relations between the Green's functions for the opposite 
directions, we can rewrite the self-consistency condition 
(J3j) separately for the real and imaginary parts of A as 

^M=7rT^Re/+K,x), (6a) 
^- = nTj2^f-(^n,x). (6b) 



Homogeneous solution for superconductor in equilibrium 
is well known, 

/| 0) = S n A RQ , /i 0) = iS n A!, = S n uj' n (7) 

with S n = l/\/< 2 + |Ao| 2 . (8) 

The absolute value of the equilibrium order parameter 
and supercurrent can be found from Eqs. ([3]) and (|4]) as 



nT y ( Re - 1 - 

tr v \/< + i A oi 2 Kiy 

j = —2eT y ] Im- 



lnT, 



(9) 
(10) 



^ 2 + |Aol 2 



Note that the depairing action of superfluid velocity in 
our model is equivalent to the depairing due to the Zee- 
man splitting, see, e.g., Ref. 1 13. As a consequence, in our 
system the temperature dependence of the critical cur- 
rent is nonmonotonic, similar to the Zeeman-effect upper 
critical field. Such temperature dependence is a peculiar 
property of clean one-dimensional superconductors and 
it is very different from higher dimensional systems. 

We will demonstrate now that Eqs. (|5|) and ([6|) allow 
for exact analytical solution. First, from Eqs. §5§ we 
obtain a useful relation between the Green's functions 



iA 



df- 
dx 



df+ , dg 



dx 



'dx 







(11) 
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In general, we cannot reduce this relation to the first in- 
tegral of equations ((5^-c). To proceed, we make some 
strong assumptions about the form of the order parame- 
ter and Green's functions. We will check consistency of 
these assumptions from the obtained solution. Namely, 
we assume that (i) the imaginary part of the order pa- 
rameter A/ is coordinate independent and (ii) /+ can be 
found in the variable-separated form 



f + (ui n ,x) = S n A R (x), 



(12) 



where the frequency term S n is assumed to be coordinate 
independent and therefore it has to be given by Eq. (|8j). 
In this case we can rewrite Eq. (fnj) as a first integral, 



d_ 

dx 



iAi/_ 



A 2 

c R 

bn ~2 



u' n g = 0, 



A 

which gives iAj/_ — S n -^ — oj' n g 
C can be found using asymptotics far away from the 
phase-slip center. Substituting for Green's functions 
their equilibrium values from Eq. (|7|), we obtain C = 



C. The constant 



-S„ [io' n 2 + Aj + A| /2) , which gives 



-Sn [ ^ + A? + *SL 



A 2 

R 



(13) 



Taking derivative of Eq. (|5c[) and excluding dg/dx and 
df-/dx from Eqs. (|5aj) and J5b|, we obtain 



d 2 f + 
dx 2 



■4 (A? 



2 )/ + -4A fl (*A 7 /_-<g) = 0. (14) 



Substituting for / + the presentation (fT2"]) and using the 
relation (|13p . we see that the frequency part indeed drops 
out and this equation reduces to the differential equation 
for A R , 



d 2 A R 
dx 2 



2A 2 , A fl 



2A\ = 0, 



(15) 



which essentially coincides with the Ginzburg-Landau 
equation. However, in our case it describes the phase-slip 
center in the whole temperature range. This equation has 
a kink soliton-likc solution 



A R (x) = A R0 tanh (A RO x 
Corresponding Green's functions are 
/+ = S n A R0 tanh (A R0 x) , 

/- = S n UAj - 



A 2 



1 



g = s„ 



2« + cosh 2 (A RO x) 

A 2 1 

^RO t 

2K + iA J )cosh 2 (A m x) 



(16) 

(17a) 
(17b) 

(17c) 



One can verify that the normalization constraint g 2 + 
/ 2 — f 2 = 1 is satisfied. The assumed coordinate inde- 
pendence of the imaginary part of A requires that the 



coordinate-dependent part of Eq. (|6b[) vanishes. This 
gives the following equation 



1 



1 



iAj 



,' 2 



= 0. (18) 



which can be used to find A/, while |Ao| is determined 
by Eq. {2} and A R0 = ^IAqI 2 - Aj. Therefore, the 
structure of the phase-slip center is fully defined. 

The probability of a phase slip event driven by thermal 
fluctuation is proportional to the exponent exp (-6F/T), 
where SF is the free-energy barrier between homogeneous 
and saddle-point state and T is the temperature. There- 
fore, the next point is to find the free-energy barrier. 

The formula expressing the free energy in terms of 
the quasiclassical Green's functions has been derived by 
Eilenberger in his original paper (l2| . This formula can 
be straightforwardly generalized to the one-dimensional 
case and in real units it can be written as 



F = N(0) / dx 



f P 



A/* + A*/ + 2u> n (g — 1) 



(19) 



The phase-slip barrier is given by the energy difference 
between homogeneous and inhomogeneous states, SF — 
F — -Fhom ; and the free-energy density of a homogeneous 
state is known to be 



^ hom = |A| 2 /A-^T 



; 2 + |A| 2 -K| . (20) 



Using the Eilenberger equations and constraint g 2 + 
fP = 1, we can simplify the expression for free en- 
ergy (fT9|) and rewrite it without space derivatives. With 
the self consistency equation we can get rid of the cou- 
pling constant A. Subtracting the homogeneous-state free 
energy from the full expression and using dimensionless 
units, we obtain the free energy in units of T c , 



SF 



5J-dx 



(21) 



with 

^ = 7r- 1 (|A| 2 -|A | 2 )lnT+ 
|A| 2 -|A | 2 _ (A 

f P 



TRe£ 



2 

S n 



(22) 



Using our saddle-point solution (TI7j|) and ((TTJ) and per- 
forming integration over x, we obtain for SF, 




A \ 2 ~iA!u;' n + A m JLj' n 2 + \Ao\ 2 



V< 2 + IA0I 2 



A RO Ju;' n 2 + \Ao\ 2 



(23) 
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FIG. 1: The temperature dependences of the phase-slip bar- 
rier in units of A(0) for different currents. The legend indi- 
cates values of supercurrent measured in units of jo = 267^. 
Note that at finite current the barrier always jumps to zero 
from a finite value. Inset shows current dependences of the 
barrier at different temperatures. 



This formula gives the barrier in terms of the sum over 
Matsubara frequencies, which allows for analytical calcu- 
lation in the limits T — > 0, 1 and can be easily computed 
numerically in general case. 

At zero temperature sum over Matsubara frequencies 
in (f2"3"| becomes an integral T^2 U ~> JdLu/2ir. At 
T — > the gap parameter |A | is given by the BCS value, 
A(0) = (tt/j)T c with 7 = 1.7811 and in our clean case 
it remains current-independent up to the critical current 
corresponding to k sa 0.83. Eq. © gives Aj — > — k/2 
for T — > 0. Integration over uo in Eq. (f2"3")l can be eval- 
uated and leads to a very simple exact result for the 
zero-temperature barrier 



We also show in the inset the current dependences of 6F 
at several temperatures. One can see that the tempera- 
ture dependences of the barrier have almost linear part 
at T — > 0. This is a consequence of the superconduct- 
ing gap vanishing at one point. It is well known that 
in fully-gapped superconductors the temperature depen- 
dences of all parameters are exponential. Accidently, the 
Ginzburg-Landau asymptotics (|25|) actually provides a 
reasonable approximation for the exact barrier in the 
whole temperature range. Even at T — > the GL re- 
sult exceeds the exact one only by 15%. 

Our results can be straightforwardly generalized for 
a case of several identical conducting channels. If we 
have M channels with the same Fermi velocity, then the 
corresponding free-energy barrier is just proportional to 
the number of channels, SFm — MSF. The current scale 
in Fig. [1] also scales with the number of channels, jo — ¥ 
]m = 2MeT c . It is clear that the thicker the wire is 
and the more conducting channels it has, the higher the 
free-energy barrier is. 

In conclusion, we considered the model of thermally 
activated phase slips in a clean one-dimensional super- 
conductor with only one conducting channel and found 
the exact analytical solution of Eilenberger equations cor- 
responding to the free-energy saddle point. This solution 
allowed us to calculate the temperature and current de- 
pendences of the free-energy barrier for thermal phase 
slips. 

The authors would like to thank I. Beloborodov for 
critical reading of the manuscript and useful comments. 
This work was supported by the U. S. DOE, Office of 
Science, under contract # DE-AC02-06CH11357. 



5F(0,fc) = -A fl0 = -A(0)l 

7T 7T 



1- 



for k < 0.83 

We also can find the asymptotics of the barrier in the 
Ginzburg-Landau region, \T— T c \ <C 1, where our results 
have to be similar to the LAMH solution. In this tem- 
perature range A<1 and we can take into account only 
the first non-vanishing term of the Taylor expansion of 
Eq. (|23p . In this limit for a wire without current we have 



SF = C A(0) (1 - T) 



3/2 



IS 



where C* = */* 1 » 0.737 with t(z) = E^-i ™~ 
the Riemann ^-function, £(3) rj 1.202. 

For arbitrary temperatures and currents the free- 
energy barrier cannot be represented in a simple form 
and has be computed numerically. Figure [1] shows the 
temperature dependences of the barrier in units of A(0) 
for different dimensionless currents in units of ja = 2eT c . 
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